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Abstract. 

We present an accurate ab initio method of calculating transition energies 
and isotope shifts in the 3rf-transition metals. It extends previous work that 
combines the configuration-interaction calculation with many-body perturbation 
theory by including the effective three-body interaction and modification of the 
energy denominator. We show that these effects are of importance in Ti II. 
The need to develop methods that can accurately calculate isotope shifts in 3d- 
transition metals comes from studies of quasar absorption spectra that seek to 
measure possible variation of the fine-structure constant a over the lifetime of the 
Universe. Isotope shift can also be used to measure isotope abundances in gas 
clouds in the early Universe, which are needed in order to test models of chemical 
evolution. 
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1. Introduction 

This work is motivated by studies of quasar absorption systems that are designed to 
probe the fine-structure constant, a — e 2 /hc, in the early Universe. By comparing 
the absorbed frequencies with terrestrial measurements one can deduce whether a has 
changed. Some studies report significant deviation from zero change (e.g. [TJ [2l [3l [4J [5] ) , 
while others do not (e.g. [6j [9] ) . 

These studies use the "many multiplet" method [10) where many transitions 
in many ions are used. This method offers an order-of-magnitude improvement 
in sensitivity over the previous "alkali-doublet" method, but introduces a potential 
systematic effect related to the isotope abundances of the absorbers: the isotope 
abundance ratios in the absorbing systems could be very different from those on earth. 
A "conspiracy" of several differences may provide an alternative explanation for the 
observed differences in spectra [TTJ [12] . 

We can resolve this problem by taking combinations of the transition frequencies 
that are insensitive to either a- variation or isotopic abundances [13) . We can then 
remove the systematic effects from the a-variation studies and simultaneously measure 
isotope abundances in the gas clouds. The measured isotope abundances can then be 
used to test models of chemical evolution in the Universe. However, to do this type 
of analysis we must know both the relativistic shift and the isotope shift of each 
transition used in the analysis. While the relativistic shifts have been calculated for 
all important ions used in quasar absorption studies, there are still some gaps in the 
isotope shift calculations because these are generally more difficult. In particular, the 
isotope shifts of the important Fe II transitions are not known. 

Previously we have calculated the isotope shift in atoms and ions with one valence 
electron using many-body perturbation theory 14j, and for many- valence-electron 
atoms and ions using a combination of configuration interaction (CI) and many-body 
perturbation theory (MBPT) [TSJ [TB]. The CI+MBPT method compares well with 
accurate multiconfiguration Hartree-Fock approaches (see, e.g. carbon calculations 
[171 fT8] compared in [IB] : Mg I calculations [19] compared in [15]), however we 
believe that our method is more readily applicable to heavier ions, hence this study 
of Ti II. This ion has additional importance because it has been observed in quasar 
absorption spectra. Combined with accurate laboratory wavelengths [20 and our 
previous calculations of the relativistic shift [21] . Ti II may be a useful probe of co- 
variation. 

In this study we extend our previous work on calculation of energy and isotope 
shift using the CI+MBPT method by including the effective three-body second-order 
MBPT operator fsection [3~T]) and varying the MBPT energy denominator fsection lBT!?]) . 
We show that these effects are of some importance in Ti II, and will probably prove 
to be of importance for all 3<i-transition metals. In section [4] we extract isotope-shift 
constants from experiment in order to benchmark our theory. This experiment is 
of particular interest due to the lack of isotope-shift data for ionised 3d elements. 
Our final results are presented in section along with predictions for astronomically 
relevant transitions. Atomic units (h = e = m e — 1) are used throughout this paper 
except where otherwise stated. 
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2. Method 



The isotope shift of an atomic transition comes from two sources: the nuclear recoil 
("mass shift"), and the finite size of the nuclear charge distribution ("field shift"). 
The mass shift is more important for light elements, while for heavy elements the 
field shift dominates. In the case of Ti II, the field shift is small (we estimate it in 
section l4~Tj) ; this paper is concerned with the mass-shift contribution, which is more 
difficult to calculate. 

Recoil of a nucleus of mass M causes a level energy shift 

\ i / i i<j 

The first term on the right hand side is known as the normal mass shift (NMS), while 
the second is the specific mass shift (SMS). We use the non-relativistic form of the 
mass-shift operator; relativistic corrections for optical transitions in light atoms are 
on the order of few percent and can be neglected [22]. We calculate the frequency 
shift of a transition between two isotopes with mass number A and A' as 



1 ' ' --'Jl\A',A 



6v A '> A = v A ' - v A 

= (%ms + ^sms) I ^7 - i ) + F6 ( r ' A ) A ' ' A • (2) 

The normal-mass-shift constant is easily calculated from the transition frequency: 
^NMS = — f/1823, where the value 1823 refers to the ratio of the atomic mass unit to 
the electron mass. The last term in this equation is the field-shift component which 
depends on the change in mean-square nuclear radius (r 2 ) and the field-shift constant 
F, calculated in section [4~T1 

To calculate the specific- mass-shift constant, £x;mSj we use ^ e a ll -orc ' er finite- 
field scaling method. Here a rescaled two-body SMS operator is added to the Coulomb 
potential everywhere that it appears in an energy calculation: 

Q = , 1 + Api • p 2 ■ (3) 

ri -r 2 



We recover the specific-mass-shift constant as 

(4) 



duo 

fcsMS = ^ 



A=0 



The operator Q has the same symmetry and structure as the Coulomb operator (see 
Appendix A in Ref. [16]). 



3. Energy Calculation 

To calculate energies, and hence transition frequencies, we use the CI+MBPT method 
[23] implemented with the atomic structure package AMBiT. This package was 
previously used to calculate isotope shifts in Mg I [T5] and carbon ions [TB]. It is 
presented in detail in [16] : here we will present only the salient points of the calculation, 
as well as some extensions that go beyond what was done in [16] (it should be noted, 
however, that these extensions have previously been proposed in [23j). 

The first step is to solve the Dirac-Fock equations for the core and valence 
electrons. From this we generate a single-particle basis set that includes the core and 
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valence states and a large number of virtual states. Then we do the full configuration- 
interaction (CI) calculation in the frozen-core approximation. Here the many-electron 
wavefunction is expressed as a linear combination of Slater determinants \I) 

where the Cj are obtained from the eigenvalue problem 

Y J H I jC, ] = EC I (5) 

j 

and H is the CI Hamiltonian. 

Core-valence correlation effects (that necessarily go beyond the frozen-core 
approximation) are included using an MBPT operator S, which is added to the CI 
Hamiltonian (see Section III of 16J ) . In this paper we calculate E to second order in 
the operator Q (Equation ||3J)), leading to the modified eigenvalue problem 

E L +E 2fflMa) c , =Kl . (6) 

./ V M M ) 

The states \M) include all Slater determinants of the single-particle basis that have 
core excitations. We further separate the MBPT operator into one-valence-electron 
and two- valence-electron parts: E^ and E^ 2 \ respectively. Goldstone diagrams and 
analytical expressions for these are given in 16J. All diagrams are included in the 
current work, including the box diagrams of E^ 2 ) which have the "wrong" multipolarity 
of the Coulomb interaction and were deliberately excluded from [16] . 

3.1. Effective three-body interaction: E*- 3 ) 

It is mentioned in 16J that there exists in the second order of perturbation theory an 
effective three-body interaction, where three valence electrons interact via the core, 
represented by a Goldstone diagram (figure [1]) with three external lines. 



Figure 1. Effective three-valence-electron interaction of E. 
a d 

—* 1 *■ 

I 

b I 



c I 



The diagrams of this type are quick to calculate: there is only one internal 
summation and no summation over virtual states. However the number of 
corresponding effective radial integrals is huge, and storage is not possible. Therefore 
we generate these diagrams as needed and seek a way of restricting them. In practice 
we only need to include orbitals from the leading configurations. For the Ti II 
transitions we're interested in this means that we only need to include E'- 3 -' in Equation 
(JHJ) when either \I) or |J) represent 3d 2 4s or 3d 2 Ap terms. In fact, the majority 
of the contribution comes from the case where both \I) and J) represent leading 
configurations. 
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3. 2. Variation of energy denominator in perturbation theory 

Another open question in the CI+MBPT method is how to deal with the energy 
denominator in Equation ©. This question is discussed in detail in [21]; see also 
[231 I25j . The basic problem is that we must make some approximation to the 
E in the energy denominator in order to generate the modified matrix elements 
before solving the eigenvalue problem and obtaining the energy spectrum E. Two 
reasonable formalisms are the Rayleigh-Schrodinger (RS) and Brillouin-Wigner (BW) 
perturbation theories. 

In [16] we used a BW method where all connected diagrams are evaluated at 
energies that correspond to the main configuration (i.e. the lowest valence energy). It 
is similar to taking an average (Dirac-Fock) energy of the states for E. In this paper 
we explore just one method of extending this, where we add a constant 5 to every 
energy denominator: 

1 1 

E-E M ~* E-Em+5 
Our idea is to make the energy E closer to the valence electron energy; to this effect 
we take 5 = E CI — E DF for the ground-state energy. In effect this "corrects" E for the 
ground state, replacing the Dirac-Fock value with the CI energy. In our case 6 rj —0.69; 
we alter it depending on the values of E CI and E DF which in turn depend on the 
particulars of the calculation, and in particular on the SMS coefficient A (Equation 
|3])). Note that this shift can be further justified on rather general grounds because 
it restores the correct asymptotic behaviour in the BW theory for a large number of 
particles [23]. 

4. Analysis of experiment 

We compare our calculations to the Doppler-free spectroscopy measurements made by 
Gianfrani et al [25] . They have made measurements of the 3d 2 As Fg /% - id 2 Ap A F°j 2 
and 3d 2 As 4 Fg/2 - 3d 2 Ap i 2 transitions for the isotope pairs 46-48 and 46-50. Their 
results are presented in table 1 of [25]; they record the following: 





Transition 


Isotope couple 


IS (MHz) 


3d 2 4s 


4 F 9/2 - M 2 Ap 4 F 9 ° /2 


46-48 


746 (8) 




46-50 


1452 (10) 


3d 2 4s 


% /2 - 3d 2 Ap *G° 11/2 


46-48 


441 (6) 




46-50 


847 (7) 



However when one looks at the raw data for these measurements (Fig. 4 in [26] ) 
it is clear that there is a misprint in the second transition: one may easily see that 
the ^G° xl j 2 isotope shift is only very slightly smaller than the 3d 2 4p 4 F°^ 2 . In fact, 
one can take a rough reading from the fitted Lorentzian curves of Fig. 4a and obtain 
£^50,46 _ 1425(150) MHz and 5v iSA6 = 725(150) MHz. 

Additionally, a mistake seems to have been made in the calculation of the normal 
mass shift. By taking the difference between the total isotope shift and residual shift in 
table 1 of [26], one obtains 236.3 MHz and 453.7 MHz for the 46-48 and 46-50 isotope 
pairs, respectively, and independent of transition energy. Our calculated values of 
NMS are shown in table [1] 
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Table 1. Extraction of the specific mass shift constant, fcg^g, from experiment. 
All isotope shift components are in MHz; fcgtvlS ' s m GHz-amu. Our calculation 
of the field-shift (FS) contribution is presented in section [4.11 





Transition 


A (A) 


Isotopes 


IS (expt.) 


NMS 


FS 


SMS 


fc SMS (GHz-amu) 


3d 2 4s 


% /2 - 3dHp 4 F° /2 


3235 


46-48 


746 (8) 


460 


43 (19) 


243 (27) 


-268 (30) 






46 - 50 


1452 (10) 


884 


103 (21) 


465 (31) 


-267 (18) 


Zd 2 As 


% /2 - MHp <G° 1/a 


3350 


46 - 48 


725 (150) 


445 


43 (19) 


237 (169) 


-262 (187) 






46 - 50 


1425 (150) 


854 


103 (21) 


468 (171) 


-269 (98) 



Due to these considerations we have used the reported values of the measured 
isotope shift for the 3<i 2 4s 4 F 9 / 2 - 3d 2 4p Fg/ 2 transition, and for the 3<i 2 4s 4 F g / 2 - 
id 2 Ap ^G"-^ transition we use our own rough reading from the data with large errors. 

In table [1] we extract fcgMS from the experimental values of isotope shift. We 
remove the NMS and field shift (FS) components (the field shift is calculated in 
section 14.11) . and then remove the mass-dependence from the residual SMS. We 
conclude from table [1] that for consistency with experiment one requires values 
of the specific-mass-shift constant of fcsMS = —267 (18) GHz-amu and —269 (98) 
GHz-amu for the 3d 2 4s % /2 - 3d 2 4p 4 F° /2 and M 2 As 4 F g/2 - M 2 4p <&° 1/2 transitions, 
respectively. 



4-1- Field shift 

The field shift component of the isotope shift is given by 

5v A '< A ' =F5{r 2 ) A '< A = F((r 2 A ,)-(r 2 A )) 

where (r 2 ) is the square-mean charge radius. We have used the rms charge radii 
for Ti II isotopes tabulated in [27], which result in 5(r 2 ) 48 ' 46 = _o. 104(45) fm 2 and 
5( r 2)50,46 = _ .251(45) fm 2 . 

We have developed two methods of calculating the field shift constant, F. The 
first is to simply vary the nuclear radius in the code (in AMBiT the nuclear charge 
has a fermi distribution) and calculate energy at each point. The field shift constant 
is then extracted as 



d{r 2 ) ' 

The second method is a scaling method where we take the difference in the 
potentials from two different nuclear charge radii and rescale it to increase the size 
of the effect: SU(r) = A ( U A (r) — U A (r)^j , where U A (r) is the nuclear potential 

of isotope A. We add 8U(r) to the original nuclear potential, perform the energy 
calculation, and extract F as 



8{r 2 ) A '' A dX ' 

This method was used in [13] for single-valence-electron ions. 

In practice, we have found that both methods give equivalent results within a 
few percent. Similarly, choice of basis sets makes little difference (our basis sets 
are discussed in the next section); any error in the FS constant is swamped by the 
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experimental error in S(r 2 ). We use F « —410 MHz/fm 2 for both the 3d 2 As 4 F 9 / 2 - 
3d 2 Ap A F°i 2 and 3d 2 4s A Fg/ 2 - 3d 2 Ap ^G^x/i transitions. The spread of our results is 
less than 10 MHz/fm 2 . 

5. Calculation and results 

The ground state of Ti II has a ls 2 2s 2 2p 6 3s 2 3p 6 3<i 2 4s configuration, and we 
are interested in transitions to 3<i 2 4p levels. A reasonable single-particle basis 
can be obtained by solving the self-consistent Dirac-Fock equations for the 
ls 2 2s 2 2p 6 3s 2 3p 6 3d 2 electrons (i.e. in the V 1 ^^ 1 approximation) and generating other 
valence and virtual levels in the potential of these electrons. 

Valence-valence correlations are included to all orders by the CI method. The 
ls 2 2s 2 2p 6 3s 2 3p 6 electrons are treated as a frozen core, and we include all single and 
double promotions from the leading configurations 3d 2 4s and 3d 2 Ap in our calculation. 
Correlations with the frozen core (including the relatively important 3s and 3p 
orbitals) are treated using many-body perturbation theory, as explained in section [3] 
Note that because the Dirac-Fock equations were not solved on the frozen core alone, 
the "subtraction diagrams" outlined in [16j must be included. There are two relevant 
basis sets: the MBPT basis that includes a very large number of virtual levels, and a 
more restricted CI basis that is a subset of the MBPT basis. 

Our calculations are performed with a relativistic single-electron basis set made 
from B-splines, similar to those developed by Johnson et al [28] . This type of basis 
was shown to be effective in calculating transition frequencies and specific mass shifts 
in Mg I and carbon ions [T51 [TB] . 

We have included in the CI all single-particle levels up to 16spdf (that is, outside 
the frozen-core we have 3s - 16s, 3p - 16p, 3d - 16d, and 4/ - 16/). This basis is 
large enough to effectively saturate the CI. For the sums over virtual levels in the 
MBPT diagrams we are able to include a much larger basis; we have used 33spdfg. 
The results are presented in table [H where we have separated the effects of the pure 
CI calculation from the effects of the MBPT operator. Note that the effective three- 
valence-electron part of the MBPT operator, E^ 3 \ has a non-negligible impact on the 
final results. This should be contrasted to other atoms (e.g. TI [23], C [16]) where they 
were omitted. In our case this is not possible because of the large overlap between the 
3s and 3p core orbitals and the 3c? valence orbitals. 

Table [3] presents a second set of results that include the modified energy 
denominator as discussed in section 13.21 As one might expect, the size of the 
MBPT contribution has been reduced because the energy denominators are now 
generally larger in magnitude. Nonetheless it is clear that the modification of energy 
denominators is a higher-order effect than the inclusion of E itself. 

In table|3]we present our calculations of isotope shift in the astronomically relevant 
transitions of Ti II, i.e. those seen in QSO spectra. The transition frequencies and SMS 
were calculated in the same manner as the transitions in table [3j with 5 = E CI — E DF . 
The difference between these results and the 5 = results was used to estimate 
accuracy. Field shift was calculated using the method of section 14.11 The results 
are presented both in MHz and kms -1 : the latter is the preferred form for use in 
astronomy (8v = SX/X x c kms -1 ). 
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Table 2. Frequencies and fcgMS of Ti II levels relative to the ground state 
(3d 2 4s 4 i*3/2 )■ Note that for this calculation we do not include modification of 
the energy denominator; i.e. <5 = 0. The last two lines show transitions where 
comparison with experiment is possible. 



Level 



CI 



£(2) s(3) Total Experiment 



3d 2 4s 4 F 9/2 
3d 2 4p 4Gg /2 
3d 2 4p *G° 1/2 
3d 2 4p 4 F° /2 



331 
28418 
28668 
29749 



lo (cm 1 ) 
-6 140 
1151 869 
1170 909 
1226 902 



-4 
273 
258 
188 



461 
30711 
31005 
32065 



393.44 
29968.30 
30240.88 
31301.01 



3d 2 4s 4 F 9/2 
3dHp 4G° /2 
3d 2 4p *G° 11/2 
3d 2 4p 4 F° /2 

4jp 4s~io 
i-9/2 - 

*9/2 - ^ 9/2 



fc SMS (GHz-amu) 



47.4 
-182.7 
-161.4 
-61.8 



-12.2 
-285.2 
-290.2 
-141.5 



1.0 
115.6 
118.4 
-131.9 



-1.5 
-22.9 
-26.9 
-10.7 



34.7 
-375.1 
-360.2 
-345.8 
-394.9 
-380.5 



-269 (98) 
-267 (18) 



Table 3. Frequencies and fcsMS °f Ti II levels relative to the ground state 
(3d 2 4s 4 F 3/2 ), including in £ our modification of the energy denominator of 
Equation J6}: <5 = E CI — E DF . The last two lines show transitions where 
comparison with experiment is possible. 

Level 



3d 2 4s 4 F 9/2 
3d 2 4p *G° /2 



3d 2 4p ^ 
3d 2 4 P 4 F° /2 



/2 



3d 2 4s 4 F 9/2 
3d 2 4p 4Gg /2 
3d 2 4p 4G° 1/a 
3d 2 4p 4 F° /2 



4 F 



9/2 



*9/2 - * 9/2 



CI 


£« 


s< 2 ) 


£(3) 


Total 


Experiment 




u> 


(cm ) 








331 


-7 


109 


-3 


430 


393.44 


28418 


837 


677 


166 


30099 


29968.30 


28668 


850 


709 


159 


30386 


30240.88 


29749 


969 


625 


109 


31451 


31301.01 




fc SMS 


(GHz-amu) 






47.4 


-9.8 


1.3 


-1.1 


37.8 




-182.7 


-185.1 


88.5 


-16.1 


-295.4 




-161.4 


-189.3 


90.5 


-18.4 


-278.6 




-61.8 


-140.0 


-29.3 


-8.4 


-239.5 












-316.4 


-269 (98) 










-277.3 


-267 (18) 



Table 4. Frequencies and total isotope shifts of Ti II transitions to the ground 
state (3d 2 4s 4 i ? 3/ / 2 ). Isotope shifts are relative to 48 Ti. 



Transition A (A) ui (cm ) 



5/ 





expt. 


theory 


(MHz) 


(km s 1 ) 


(MHz) 


(km s 1 ) 


3385 


29544 


29636 


-784 (79) 


0.265 (27) 


742 (73) 


-0.251 (25) 


3243 


30837 


30937 


-756 (98) 


0.245 (32) 


717 (90) 


-0.232 (29) 


3230 


30959 


31075 


-742 (95) 


0.240 (31) 


704 (88) 


-0.228 (28) 


3074 


32532 


32732 


-744 (107) 


0.229 (33) 


706 (99) 


-0.217 (30) 


3067 


32603 


32821 


-743 (113) 


0.228 (35) 


705 (104) 


-0.216 (32) 



3d 2 4p 4G| /2 
3d 2 4p 4 F° /2 
3d 2 4p 4 F° /2 
3d 2 4p 4 D° /2 
3d 2 4p 4 D° /2 
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6. Conclusion 

In this paper we have calculated isotope shifts for Ti II as a test case for 3d-transition 
metals because it has experimental data available for comparison. Our CI calculations 
for Ti II designate 3c? as a valence orbital while keeping the 3s and 3p orbitals in the 
frozen core. This provides a great saving for the size of the CI calculation, but the 3s 
and 3p correlations are large, and must be included using MBPT in order to obtain 
good accuracy for the specific-mass shift. In this case we have shown that S^ 3 ^ is 
important; we would expect it to be for all 3e?-transition metals. 

We have also presented a modification of the MBPT energy denominator which 
consists of adding a constant, 5 — E CI — E DF , to the denominator. Here E CI 
and E DF are calculated for the ground state. This improves both the calculated 
frequencies and specific-mass-shift constants significantly. We have performed isotope- 
shift calculations for astronomically important transitions that are seen in quasar 
absorption spectra. Taken together with existing accurate laboratory wavelengths 
and relativistic-shift calculations, Ti II may be a useful probe of a- variation. 

Using the techniques presented in this paper we can perform accurate calculations 
of isotope shifts for other atoms of astronomical interest, including the 3d-transition 
metals Fe II, Cr II, Ni II, and Mn II. These systems are very important in quasar 
absorption studies of a-variation, yet their isotope shifts have not been measured. 
Furthermore, by comparing results of different methods, in particular the 5 = and 
S = E CI — E DF results, we can estimate the accuracy of our theoretical predictions. 
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